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SYSTEM AND METHOD FOR MODELING BONE STRUCTURE 

PRIORITY INFORMATION 

This application claims priority under 35 U.S.C. § 1 19 from US Provisional Patent 
Application Serial No. 60/240,884, Hied October 17, 2000, and US Provisional Patent 
Application Serial No. 60/246,198, Hied November 6, 2000; each of which is hereby 
incorporated by reference in its entirety. 

FIELD OF THE INVENTION 

The present invention discloses a structural and mechanical model and modeling methods 
for human bone based on bone's hierarchical structure and on its hierarchical mechanical 
behavior. The model allows for the assessment of bone deformations and the computation of 
strains and stresses due to the specific forces acting on bone during function. The model further 
contemplates forces that do or do not cause viscous effects and forces that cause either elastic or 
plastic bone deformations characterized by fractures. In preferred embodiments the model is 
computerized, for example using computer simulation, imaging and rendering techniques. 

BACKGROUND OF THE INVENTION 

Material science defines the structural properties of an object as the properties that 
describe the object's makeup independent from its shape. Adult human bone has a complex 
structure and can be described as a four order hierarchy, arranged in decreasing size (Petersen, 
1 930). The first order, macrostructure (Fig. 1), comprises the structures corresponding to gross 
shape and differentiation between compact (or cortical) bone (Fig. 2) and spongy (cancellous or 
.trabecular) bone (Fig. 3). Compact bone is present in the long bone shaft (or diaphysis). Spongy 
bone is present in the lower jaw (mandible), in the epiphysis of long bone shaft, and in flat and 
short bones. The second order (or microstriicture) of compact bone includes lamellar systems 
(lamellae). Organized lamellae around vascular canals are referred to as osteons (harvesian 
systems) and disorganized lamellae among osteons are referred to as the interstitial bone. The 
second order also is comprised of related structures such as bone marrow (see e.g. Bloom and 
Fawcetts, 1986). The third order (or ultrastruclure) ofcompact bone consists mainly of collagen 
bundles and hydroxyapatite crystallites; mucopolysaccarides amount to a small amout but may--. 
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have a significant role. The fourth order of compact bone consists of molecular arrangements 
between organic and inorganic substances. For cancellous bone, the second order includes 
trabecular which comprise lamellar systems and related structures, e.g. bone marrow. The third 
and forth orders of cancellous bone are the same those described for compact bone. 

The osteon comprises a haversian canal with concentrically arranged lamellae. Osteons 
of long bone are generally directed along the long bone axis. Osteonic lamellae are organized as 
consisting of an organic framework (mostly a collagen bundle) embedded in ground substances, 
such as proteins and water, and hydroxyapatite crystallites. The hydroxyapatite crystallites are 
oriented in directions analogous to those of the bundles. Osteons measure a few centimeters in 
length and are between 200 and 300 pm in diameter. The degree of osteon calcification (relative 
amount of hydroxyapatite crystallites) is variable from osteon to osteon as well as within osteons. 
These differences are proposed to be due to the process of bone renewal or remodeling. In this 
process, osteons are renewed continuously. Consequently, osteons at different degrees of 
calcification are always present in adult compact bone. 

There is a spectrum of osteon types that refer to the arrangements of fiber bundle direction 
in the lamellae. Two osteon types, "longitudinal" and "alternate 11 , are representative of the two 
ends of the spectrum. Longitudinal osteons consist of bundles with a marked longitudinal spiral 
course. Alternate osteons consist of bundles with a marked longitudinal, oblique, and transversal 
course in successive lamellae (Frasca et al y 1 977; Giraud-Guille, 1 988; Ascenzi A. et aL, 2000). 
There are two types of lamellae, termed extinct (or longitudinal) and bright (or transverse or 
circulariy-fibered) lamellae. Extinct (or dark) lamellae appear extinct (or dark) whereas bright 
lamellae appear bright under a polarizing microscope when the microscope and osteon axes are 
aligned. 

Compact bone consists of about 40% minerals, 40% collagen, and 20% fluids. The major 
internal spaces or discontinuities of compact bone include the vascular system, pits and cavities 
(lacunae), narrow channels (canaliculae), fine porosity, and spaces between the mineral phases. 
The major internal material discontinuities of compact bone (Fig. 5), in order of decreasing size, 
are: 



Vascular system 


20-50 fim 


Lacunae 


4-6 \x\n 
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Canaliculae 


0.5-2 (.im 


Fine porosity 


600-800 A 


Spaces between mineral phases 


50-100 A 



Cancellous bone consists of trabecular i.e. osseous structures with either a sheet-like or a 
rod-like configuration. These structures interlace to form a lattice-like or spongy biological 
structure (Fig. 3). For example, both types of trabeculae are present in the calcaneous; however, 
up to 3% of the rod-like configurations are tubular due to the vascular canal running through 
them. Therefore, they are similar to the harvesian system. In general, tubular trabeculae appear 
to have a relatively simple structure. Collagen fibrils run mostly parallel to the long axis of 
tubular trabeculae in the trabeculae outer portion and perpendicular in the inner portion. 
Although the true density of fully calcified cancellous bone is a little lower and the proteoglycan 
content a little greater than those of the fully calcified compact bone, the substantial difference 
between compact and cancellous bone resides in the porosity. The cancellous bone porosity, 
which ranges from 30% to more than 90%, is mainly due to the wide vascular and bone marrow 
intrabecular spaces. As is seen in compact bone, levels of calcification vary from trabecula to 
trabecula and within trabeculae. 

The connections and orientations of trabeculae are found to have precise patterns, which 
are believed to relate to specific mechanical properties. The structure of the cancellous bone in 
the head and in the neck of the femur is usually given as an example of the correlation between 
the orientation of the trabeculae and the linear distribution of the principal forces during load 
bearing (stress trajectoral theory (Bell, 1956)). In general, such correlation between the 
orientation of the trabeculae and the linear distribution of the principal forces during load bearing 
is still understudy because while in line with the mathematical calculations, the possible effect of 
muscle traction is complex (Koch, 1917; Rybicki et aL s 1972). Nevertheless, there is a close 
relationship between the number and arrangement of trabeculae and the strength of cancellous 
bone (see e.g. Kleerekopere/a/. ? 1985). This is evidenced by the age-induced loss of trabeculae 
(see e.g. Birkenhager-Frenkel et al, 1988). Since this loss is rather selective {i.e. transverse 
trabeculae disappear more frequently than vertical ones in the central zone of the osteoporotic 
vertebral body; entire trabeculae totally disappear in elderly women and a sharp fall in trabecular 
number is observed in elderly men), it is possible that cancellous bone contains some bundles of 
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trabeculae whose main function is io resist mechanical forces while others have mainly a 
metabolic role. 

The mechanical behavior of an object, or the response of an object to forces, of an object 
depends on the structure of the object. If the object is comprised of a hierarchical structure, the 
mechanical behavior of the object varies from order to order. That is, each order or level of the 
hierarchy responds to forces according to the structures and relationships within that order. 
Overall mechanical behavior of the object is ultimately determined by the mechanical properties 
of the different orders. Therefore, the mechanical properties of an object will vary with the 
hierarchical structure of the object. Bone is an example of an object where the mechanical 
behavior and mechanical properties are dependent upon this kind of hierarchical structure. 

Mechanical properties of bone have been and are being investigated at various 
hierarchical levels through invasive (sample isolation) and non-invasive testing. Osteonic 
trabecular lamellae, osteons, trabeculae, and macroscopic compact and cancellous bone samples 
have been and are the objects of such studies. Micromechanical results include Ascenzi A. and 
Bonucci, 1964, 1967; Ascenzi A. and Bonucci, 1968, 1972; Currey, 1969; Ascenzi A. et at., 
1985, 1997, 1998;HohIing£/aA, 1990; Ascenzi A. et at., 1990, 1994; Marotti elaL, l994;Zive/ 
al. 9 1996; Ascenzi M.-G., 1999a, 1999b; Hujaefa/., 1999;Zyssete/a/., 1999; Ascenzi M.-G. et 
aL, 2000. Macromechanical results include Hazama, 1956; Cook and Gordon, 1964; Carter and 
Hayes, 1976 and 1977; Carter et aL, 1976 and 1981; Carter and Spengler, 1978; Hayes and 
Carter, 1979; Burr et aL, 1988; Cater and Carter 1989; Jepsen and Davy, 1997. 

Even though numerous publications have addressed bone micromechanics in recent years, 
many biomechanical issues relating to bone are still not understood due to the lack of reliable or 
predictive models. The lack of inclusion of such micromechanical properties in current models 
of bone functions and behavior have severely limited their usefulness in predicting 
macromechanical properties. These properties include the bone behavior in response to external 
forces or identifying the requirements of bone reconstruction and prosthesis. However, the 
inclusion of these factors requires the development of methods and studies that may provide 
reliable and reproducible results. 

The present invention describes a method to understand and predict the behavior of bone. 
The method includes a model of macroscopic bone which is constructed in terms of bone's 
hierarchical structural and mechanical properties and their interaction with forces acting on the 
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macroscopic bone, including forces associated with the ordinary functioning of the body and 
forces applied clinically. The method can be applied to any bone structures, including human 
bone and the bones of vertebrates in general. The model applies to normal bone, and to 
pathological bone, when the pathology either does not alter the structural hierarchy, or when the 
5 alterations are characterized. The model is also applicable to fossilized bone. 

SUMMARY OF THE INVENTION 

The present invention contemplates a model of macrostruotural properties of bone. The 
model comprises hierarchical structural and hierarchical mechanical properties of microstructure 
1 0 of the bone and includes interactions of the bone with internal and external forces. In a preferred 
embodiment, the bone that is modeled is either compact bone or cancellous bone. In an 
additional preferred embodiment, the mechanical properties used in the model are selected from 
the group consisting of tension, compression, shear, bending, torsion, prestress, pinching, and 
cement line slippage. 

15 The present invention also contemplates methods of predicting deformation and 

fractures of bone and for identifying the requirements of bone reconstruction and prosthesis using 
the model of the present invention. 

BRIEF DESCRIPTION OF DRAWINGS 
20 Figure 1. A schematic representation of the upper third of the tibia; i.c.s. and o.c.s. stand 

for inner and outer circumferential systems, respectively. Both compact and cancellous bone are 
represented. 

Figure 2 (a) and (b). (a) Diagram of a diaphysis sector of cortical long bone. The 
osteons or haversian system (HA) are located between the outer OL and inner IL circumferential 
25 lamellae. Theosteonic lamellae are disposed cylindrically around the haversian canal (HC). (b) 
Cross-sectioned osteons as seen (A) under a light microscope; (B) in a microradiography and (C) 
under the polarizing microscope. 

Figure 3 (a) and (b). (a) Section of the body of a lumbar vertebra showing vertical and 
horizontal trabecular The upper and lower surfaces correspond to articular cartilage, (b) 
30 Section of half of tibia's upper third. The cancellous bone of the metaphysis consists of 
comparatively think vertical trabeculae connected by thin trabecular 
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Figure 4 (a)-(c). (a) Cross section of an isolated longitudinal osteon, magnified 270 
times, (b) Cross section of an isolated alternate osteon, magnified 270 times, (c) An isolated 
osteonic sample with lugs, magnified 20 times. Lugs are used to grab the sample during 
mechanical testing. Dimensions: inner diameter 52 pin, outer diameter 225 j.un, length 500 yon. 

Figure 5. Cross-section diagram of an osteon sample illustrating the arrangement of 
canaliculae and lacunae relative to lamellae. 

Figure 6 (a)-(c). (a) Types of pure forces, (b) Definition of stress on an area on which the 
force is constant, (c) Definition of unidirectional strain for D much smaller than L. 

Figure 7(a) and (b). (a) Tensile and compressive stress distribution during torsion in a 
material, such as macroscopic bone, which is weaker in tension than in shear, (b) Shearing stress 
on the cross section of a specimen subjected to torsion. The arrows' length indicates the 
magnitude of the shearing stress, which progressively increases from the center to the periphery 
of the specimen. 

Figure 8 (a) and (b). (a) Bending of femur due to gravity, C indicates the area under 
compression and T indicates the area under tension, (b) Diagrams A, B, and C display the 
distribution of transverse and longitudinal lamellae in the sections prepared from the upper, 
middle and lower shaft, respectively. The posterior, anterior, medial and lateral regions 
correspond to the top, bottom, left and right regions, respectively, on the page. The distance 
between the centers of two adjacent square symbols measures 1 .86 mm. The size of the square 
symbol is proportional to the ratio of the bright area in circularly polarized light to bright area in a 
dark field illumination. The regions with dominant transverse lamellae correspond to the regions 
with concentration of larger squares in the upper medial, middle medial-posterior and lower 
posterior shaft, which correspond to the areas of compression in (a). The regions with dominant 
longitudinal lamellae correspond to the regions with concentration of smaller squares in the upper 
lateral, middle lateral-anterior, and lower anterior shaft, which correspond to the areas of tension 
in (a). 

Figure 9 (a)-(c). (a) The osteonic lamellar model is a laminate, which consists of fiber- 
reinforced unidirectional laminae, (b) The interstitial lamellar model is a portion of the osteonic 
lamellar model. The figure shows three thin laminae (lamellae) and a thick lamina (portion of 
cement line), (c) On a small laminar element of constant thickness, the principal material axes are 
labeled 1, 2, and 3. Direction 1 is parallel and direction 2 is perpendicular to the fibers. 
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Direction 3 is the radial direction perpendicular to the page. Circumferential and axial directions 
are labeled 0 and z. The angle between the circumferential direction and direction 1 is called y. 

Figure 10. Shows a device for subjecting bone to torsional cyclical loading. 

Figure 11. A schematic diagram of a device for subjecting bone to torsional cyclical 
loading, where (1) is a rotational axis with jaws; (2) and (3) are hard metal wedges of a pendulum 
loading system; (4) is a wheel around which a tungsten thread loaded with weights is attached; 
(5) is the axis of the pendulum; and (6) is a mirror. 

Figure 12. An idealized bilinear hysteresis model of curve prior to cycling and a first 
cycling loop, where (a) pinching is present; and (b) pinching is not present. 

Figure 13. A diagram that shows that around each osteon sample, a trapezoid was cut 
with a blade under a stereo microscope. 

Figure 14. A diagram that shows that after isolation, each lamellar sample was carefully 
straightened to a ribbon-like shape. 

F12 ure 15. A larger view of the lamella described in Figure 14. 

Figure 16. A lamella after tensional testing. 

Figure 17. Collagen bundles of bright lamella under polarized microscope. 

Figure 18. Collagen bundles of bright lamella under polarized microscope. 

Figure 19. Collagen bundles of extinct lamella under polarized microscope. Bundles are 
parallel to the osteon axis when embedded in bone. 

Figure 20 (a)-(h). Isolated and flattened bright lamella under the confocal microscope. 
From border to center, the collagen bundles go from oblique to vertical. 

Figure 21 (a)-(g). Isolated and flattened extinct lamella under the confocal microscope. 
From one border to the other, the collagen bundles are parallel to the osteonal axis. 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention describes a method for modeling the anisotropic (direction- 
dependent) and non-homogeneous macrostructural properties of compact bone in terms of the 
microstructure. The model is based on the hierarchical structural and mechanical properties and 
bone interactions with internal and external forces. An example of such forces includes, but is 
not* limited to, the ordinary functioning of the body. The model includes properties of the 
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microstructure, in particular distributions of transverse lamellae of trabeculae and of alternate 
osteons, pinching of osteons, and slippage of osteons at the cement line. 

Morphological and mechanical studies of bone show that at all hierarchical levels bone is 
anisotropic (the local mechanical properties are direction dependent), and non-homogeneous (the 
structure is not the same at different points). Nevertheless to simplify bone modeling, bone 
structure often is assumed to be homogeneous, isotropic (not direction dependent), transversely 
isotropic (one plane of symmetry), or orthotopic (three planes of symmetry). The simplifications 
of isotropy, orthotropy, and transverse isotropy give rise to unrealistic models because these 
simplifications that symmetries that do not exist. For instance, in such models stresses may be 
over- or under-estimated. When such models are applied to practical applications, for example 
bone implants, poor estimates of stress may give rise to screw loosening in implants. The 
simplification of homogeneity gives rise to unrealistic models because it disregards the hierarchy 
of the bone structure. The existing hierarchical models are based on homogeneity theory, finite 
element analysis, and classic and Cosserat elasticity theories (see e.g. Kate, and Meunier, 1987; 
Crolete/a/., 1993; Pidaparti and Burr, 1992). These models do not include important properties 
of the microstructure, which are included in the present invention and are described below. 

The present model provides for modeling each level of the hierarchical structure of bone 
in terms of the structural and mechanical properties of that level. The model further provides for 
determining the relationships among the various levels. 

In compact bone: 

(1) Collagen bundles, hydroxyapatite -crystallites, and mucopolysaccharides are 
organized in two lamellar types, bright, which are prestressed, and extinct lamellae. 
Lamellae show porosity. 

(2) Lamellae are organized in osteons. Osteons show pinching under tension- 
compression cyclic loading. 

(3) Osteons samples are organized in osteon sample groups on bone sections. 
Osteon groups show cement line slippage during torsion. 

(4) Osteon groups are organized to complete a bone section. The collagen bundle 
direction distributions are used to complete this organization. 

(5) Transverse sections are organized to complete the macroscopic bone. 
In cancellous bone: 
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(1) Collagen bundles, hydroxyapatite crystallites, and mucopolysaccharides are 
arranged in bright and extinct lamellae. Bright lamellae are prestressed. Lamellae show 
porosity. 

(2) Lamellae fill in trabecular Trabeculae show porosity. 

5 (3) Trabeculae are grouped in trabecular sample groups on bone sections. 

(4) Trabecular groups are organized to complete a bone section. The prevalent 
collagen bundle directions are used to complete this organization. 

(5) Sections are organized to complete the macroscopic cancellous bone. 

10 To produce the model, each microstructural level of the hierarchical bone structure 

commences with the micro-structural components and proceeds through the macro-structure. 
Each element of the assembly is correlated with mechanical properties either determined from 
literature sources or that are newly estimated. Homogenization methods are used to assemble the 
structure at one level with the structure at the next level, and so on, to build a hierarchical model. 

1 5 The finite element method allows for the computation of strains and stresses throughout the 
model. 

The invention results in a model of macroscopic compact and cancellous bone that 
respects the hierarchical, structural, and mechanical properties starting from the microstructural 
components. The model may be applied to all bones to result in a model for each bone of the 
20 skeleton. 

The present invention further defines methods of predicting deformation and fractures of 
bone and identifying the requirements of bone reconstruction and prosthesis using the model. 
From the specific forces that act on bone during function, the model allows for the assessment of 
bone deformation, strain, stresses, and fractures. Additionally, from the fractures and stress 

25 distribution, the model allows for the computation of strain deformation and forces that cause the 
observed fracture and stress distributions. The model also contemplates forces that do or do not 
cause viscous effects. The model contemplates forces that cause either elastic or plastic bone 
deformations as characterized by fractures. 

The model includes torsional cyclic loading functions of two representative osteon types 

30 in terms of degrading properties such as stiffness and pinching, and increasing energy absorption. 
These mechanical property changes are correlated to the idealized or mathematical behavior of 
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ultrastructural components, which includes yielding, buckling, and fracturing properties. The 
resulting algorithms and behaviors comprise an osteon model, which simulates fracture 
propagation in osteons under cycl ic torsional loading in terms of microcracking, debonding, void 
growth, and fiber breakage. Verification of the model is demonstrated by checking that the 
model produces the fractures observed in osteon samples that are separately submitted to tension, 
compression, and shear. 

The model also includes simulation of micros true tural fracture propagation in bone. 
Because of the dependence of the macrostructure's mechanical properties on microstructure, the 
model will provide an improved understanding of properties of long bone, such as fracture 
propagation, including a better understanding of how human bone macrostructure responds to 
forces acting on it. 

The model will have application in many areas, including without limitation: 

• the mechanics of natural composites and the manufacture of new composites, since bone is a 
natural composite material; 

• the identification of the fundamental requirements of bone reconstruction and prostheses 
(which will increase design effectiveness and reduce testing and related cost); and 

• the microstructure of vertebrates whose microstructure is similar to human's 

Definitions 

The present invention spans through both elastic and plastic ranges. As used herein, the 
term "elastic range" refers to the stress and strain values for which the material structure does not 
break and returns to its original shape when the force is removed. As used herein, the term 
"plastic range" refers to the stress and strain values for which the material structure does break 
and therefore does not return to its original shape when the force is removed. When an 
increasing force (starting from zero) is applied to a material, the material undergoes first elastic 
and then plastic deformation. Any bone type can undergo elastic deformation only or both elastic 
and plastic deformation depending on the force magnitude. Elastic and plastic deformations 
provide a starting point to predict strain and stress distributions and fractures of bone. The model 
also may be used to compute the stress distribution from the strain distribution and strain 
distribution from elastic and plastic deformations. It" further identifies the requirements of bone 
reconstruction and prostheses. 
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As used herein, the term "boundary conditions" refers to the relative movements of the 
boundaries of the various hierarchical structures under loading. In a specific embodiment, the 
behavior of the cement line under loading is the boundary conditions for the osteon and the 
interstitial bone between which the cement line lies. 

The term "pinching" refers to a sharp change of stiffness of bone. As used herein, the 
change in stiffness can be either from increasing stiffness to decreasing stiffness or from 
decreasing stiffness to increasing stiffness. In a preferred embodiment, the change presents itself 
on each half-cycle. 

The term "material analog" refers to a model or reproduction produced from material, as 
distinguished from a mathematical or computer model. 

The term "distraction device" refers to an apparatus that generates bone by stimulating 
growth of existing bone by application of forces to such existing bone. 

As used herein, the term "strain distribution" refers to a measure of the degree of 
elongation at any point on a sample. In a preferred embodiment, the sample is bone. 

The term "stress" refers to the force per unit area. 

The term "stress distribution" refers to strain distribution and on the mechanical property 
distribution throughout the body. 

The term "corrected break area" refers to the actual bone area, except for the lacunae and 
canaliculi, subjected to stress in the vicinity of a break. 

The terms "viscous effect" and "viscoelastic" refer to a system that exhibits behavior that 
combines liquid-like and solid-like characteristics. 

Factors 

Mechanical Properties 

Various mechanical properties are included in the model of the present invention. The 
properties will be correlated with each hierarchical level of the bone to produce the model. A 
non-limiting list of such properties is disclosed and described herein. 

The mechanical properties of bone are quantified by parameters or coefficients that 
describe the response of bone to tension, compression, shear, bending and torsion. Tension, 
compression, and shear are termed "pure forces" because each of them is recognized by the effect 
(the deformation) produced in the body to which it is applied (Fig. 6). A tension (tensile) force 
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tends to lengthen the body to which it is applied, while a compression (compressive) force has a 
tendency to shorten the body. A shear force tends to make one part of a body slide in a direction 
opposite to that of an adjacent part. Bending and torsion (Fig. 7) are a combination of tension, 
compression and shear. 

The effect of the application of one of the above-mentioned forces to a body at a natural 
state is described in terms of strain and stress. Strain is the measure of dimensional changes in a 
body and is computed by means of the deformation (Antman, 1995). Since in general the value 
of strain changes from point to point throughout the body, more properly one refers to strain as 
the strain distribution throughout the body, which provides the value of strain at each pointofthe 
body. The tendency of a body to be deformed by the application of a force is resisted by the 
internal force among the molecules composing the body. Such resistance is measured by the 
stress, which is a force per unit area. Similar to strain, in general the value of stress changes from 
point to point throughout the body. More properly one refers to stress as the stress distribution 
throughout the body, which provides the value of stress at each point of the body. The stress 
distribution depends on the strain distribution and on the mechanical property distribution 
throughout the body. What all elastic structures have in common is that the stress distribution is 
a linear function of strain within the elastic range (Hooke's Law; see e.g. Jones, 1975). Beyond 
the elastic range the relationship between stress and strain distributions depends on the particular 
structures. For instance, the mechanical testing of a specimen provides a stress-strain diagram, 
which allows the study of the relationship between stress and strain. 

Studies indicate that the mechanical behavior of longitudinal and alternate osteon samples 
at equal degree of calcification, as assessed by the method of Amprino and Engstrom (1952), 
differs because of their structural difference. The comparison of experimental stress-strain 
diagrams for longitudinal and alternate osteons shows that under monotonic tension and torsion, 
longitudinal osteon samples resist stresses better than alternate osteons; while under compression, 
shearing and bending, alternate osteon samples resist stresses better than longitudinal osteons. 
Under cyclic tension-compression loading, longitudinal osteon samples show a larger energy loss 
and lesser pinching degradation per cycle than alternate osteons; longitudinal osteon samples 
show a greater strain increase during compression than tension. The opposite is true for alternate 
osteons. ' . 
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The macroscopic mechanical properties have been found to depend on and to be 
explained by the microstructure. In particular, they have been found to depend on the numerical 
presence of osteons, the size and percentage volume of osteons, and collagen fiber orientation 
(Currey, 1959; Evans and Vincentelli, 1969). As early as 1873, Rauber considered the 
correspondence between bone micro- and macro- structure. He hypothesized that the structure of 
osteons and interstitial bone in the long bone shaft relates both to their distribution in the shaft 
under normal conditions, and also under pathological conditions that do not alter the bone's 
hierarchical configuration. This hypothesis was later confirmed (Portigliarti-Barbos, 1983, 1984, 
and 1987; BoydeeM/., 1984; Ascenzi A. eta/., 1 987a and 1987b; Ascenzi A., 1988; Carandoe/ 
a!., 1989 and 1991). Specifically, the distribution of dark lamellae (whose bundles have a 
transverse and oblique course) and of bright lamellae (whose bundles have a longitudinal course) 
in osteonic and interstitial bone follows a characteristic non-random pattern. Studies indicate that 
this distribution is consistent with the distribution of bending forces usually operative on this 
bone (Ascenzi M.-G., 1999a). For example, the femoral dominant distribution of dark and bright 
lamellae displays a clockwise rotation of approximately 90° in sequential sections from upper, 
middle, and lower third of the shaft (Portigliatti-Barbos, 1983, 1984). In fact, because of the 
femoral overall shape which includes two curvatures (an anterior-posterior curvature and a 
lateral-medial curvature), gravity on the body results in the bending of the femur. Because 
bending always includes an area of tension and an area in compression, the femur presents an 
area in tension and an area in compression (Fig. 8a). It turns out that the femoral dominant 
distribution of dark (bright, respectively) lamellae coincides with the area in tension 
(compression, respectively) (Fig. 8b). Recent work also has found that the above-mentioned 
transverse/longitudinal lamellar distribution is consistent with the distribution of alternate osteons 
(Herte/a/., 1994). Neither the transverse/longitudinal lamellar distribution nor the alternate 
osteon distribution have been included in prior models of bone structure. Transverse/longitudinal 
lamellar distribution and the alternate osteon distribution are included in exemplary models of 
the present invention. 

Additional Factors 

The following provides a non-limiting list of factors that may be included in the models 
of the present invention, and which are used in exemplary embodiments. 
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(I) Fracture of macroscopic bone. The invention incorporates fracture dynamics into the 
bone model and modeling methods, including mechanisms by which a fracture starts and spreads. 
Unlike other models of bone, fracture propagation is modeled in terms of ultrastructural 
components. The literature indicates that the fracture mechanism of bone depends on bone 
structural and composition properties such as collagen architecture and collagen content (e.g., 
Jepsen ei a/., 1999). In 1969, Evans and Vincentelli showed significant differences among 
osteons of various bones (fibula, tibia, and femur) in the "corrected break area", which is the 
actual bone area, except for the lacunae and canaiiculi, subjected to stress in the vicinity of a 
break. Characteristic differences were found between the means of the corrected break area for 
groups of longitudinal and transverse {i.e., consisting of transverse lamellae) osteons and of 
osteon fragments of the femoral and tibia sections and for groups of the transverse osteons and 
fragments of the tibia and fibula sections. The percentage of the "corrected break area'* of 
transverse osteons and their fragments in the tibia and fibula sections was also statistically 
different. Another study (Vincentelli and Evans, 197 1) established a relationship among macro- 
mechanical properties, collagen bundles, and calcification in the shaft of long bones. 
Furthermore, fracture lines appear to follow the cement lines between osteons and lamellar 
boundaries within osteons (Simkin and Robin, 1 974) where the bone is weaker. According to the 
invention, inclusion of the differences between the means of the corrected break area for groups 
of osteons would increase the predictability of the present model compared to prior art models. 

(2) The prestress distribution in bone. The models and methods of the invention 
incorporate computations of the stress distribution in long bone so as to include pre-stress 
(Gurrey 1964; Ascenzi A. and Benvenuti, 1980). Stress distribution in long bone depends on 
structural and composition properties such as collagen architecture and collagen content. Bone 
areas where collagen bundles are transverse and oblique to the long bone axis are prestressed. 
Such prestress, estimated on the order of 0. 1 GPa, is too large to be disregarded. It locally 
impacts the stress produced by forces acting on bone (see Ascenzi M.-G., 1999a). Newly 
estimated prestress variables are included in the present invention (Ascenzi M.-G., 1998a and 
1999b). The newly estimated prestress was evaluated through the structural and mechanical 
modeling of isolated lamellar samples, and has been shown to be a realistic approach. See e.g., 
A. Meunier, 1 999. Inclusion of this prestress into the model of the present invention allows one 
to more accurately model bone in terms of computation of stresses. 
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(3) The phenomenon of "pinching 11 . The invention for the first time incorporates 
pinching into bone models and modeling methods. Pinching is the mechanism of yielding and 
buckling of collagen bundles under loading beyond the elastic phase. It is an important step in 
the formation and propagation of fractures. The understanding of pinching requires detailed 
analysis of osteon mechanical behavior. In fact; while the stress-strain curves for monotonic 
loading under tension, compression, and torsion show trends no different from those recorded 
from macroscopic bone samples, the tension-compression hysteretic loops showed a new 
behavior for osteon samples not observed in macroscopic samples (Ascenzi et al, 1985 and 
1997). The new behavior observed is that tension-compression hysteretic loops of osteons 
demonstrate S-shaped half-cycles. This phenomenon has been observed and studied only relative 
to earthquake-resisting structures. In such context the behavior is usually called "pinching" (see, 
e.g. Narayanan and Roberts, 1991). Pinched hysteretic loops are typical of structures that 
incorporate a matrix that cracks and reinforcements that yield or whose members buckle when 
subjected to compressive loads. In osteons, the shape and dimensions of hydroxyapatite 
crystallites and the relationship of these parameters to the organic components of the matrix are 
only partially known. Not all collagen bundles are completely calcified. Those that are not 
calcified take up crystallites only on 400 A bands (Ascenzi A. et al, 1 965). Hence such bundles 
may be comprised of relatively more stiff 400 A bands separated by relatively more flexible non- 
calcified collagen segments. 

Pinching in osteons is hypothesized to be mainly localized at the partially calcified 
bundles. Therefore, in osteons, either bundles yield in tension and buckle in compression while 
crystallites fracture and detach from collagen, or crystallites fracture and detach from collagen in 
tension while collagen yields in compression. Thus, cyclic tension-compression loading shows 
pinching. Since cyclic torsional loading involves tension and compression, cyclic torsional 
loading is expected to show pinching. Nevertheless, it may be that the disruption created by 
torsional loading is too disordered, in comparison to that due to tension-compression, to allow 
closing of lesions and resolution of members as controlled as under tension-compression. In any 
event, if cyclic torsional loading of osteons shows pinching, pinching is included in the invention 
applied to macroscopic compact bone torsional loading. 

(4) Macrosintctiire and mechanical loadings/tidies of whole bone or macroscwiples. The 
invention takes account the influence of bone microstructure in evaluating mechanical loading of 
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whole bone and of bone macrosamples. In the literature, for example the torsional loading in 
bone has been analyzed using finite element analyses (see e.g. Hazama, 1956; Pfafrod etal, 1972 
and 1975; Kneise/aL, 1973; Miller and Piotrowski, 1974; Evans, 1 978; Martens ei al, 1980; 
Moreland, 1 980). However, models presently do not completely reflect the changing properties 
of bone at the mici ostructural level. Similarly, cancellous bone has been described as continuous 
and isotropic, which does not reflect the high porosity and the changing details (such as collagen 
bundles direction and lamellar structure) at the microstructural level. The elastic and plastic 
moduli change locally in relation to the microstructural properties. 

Such studies ignore most of the mechanical properties of the microstructure (because 
macroscopic samples do not always have the same mechanical properties as the microstructure 
that comprises them) and therefore do not provide a realistic understanding of bone mechanics. 
For instance, pinching is present in longitudinal and alternate osteons but not in macroscopic 
compact bone samples during tension-compression cyclic loading; also the torsional stiffness 
varies from osteon samples to osteon groups and relative to that of larger compact bone samples 
(Lakes, 1995). Lakes shows that the torsional shear moduli of osteons are much larger than shear 
moduli obtained for macroscopic samples. That is, slender specimens are stiffer than thick ones; 
the lower stiffness in thick specimens is attributed to slippage of osteons at the cement lines 
during torsion of macrosamples. Such slippage is described well by Cosserat elasticity theory 
since it allows a moment per unit area in addition to the usual force per unit area of classic 
elasticity theory. The inclusion of this factor into the model impacts, for example, the simulation 
of fracture propagation. The fracture propagation model is able to simulate the slippage of 
osteons at the cement lines during torsion and therefore the experimentally obtained results 
regarding fractures spreading along the cement line. 

Local Properties and Bone Modeling 

The knowledge of the mechanical properties and of the strain and stress distributions of 
compact and cancellous bone under specific loading is necessary in all contexts where the local 
behavior of bone is in question. For example, stability is the crucial characteristic of an 
osteotomy fixation device. When a tibia requires an osteotomy, the device that holds in place the 
two bone edges created by the cut (osteotomy) can only alfow for micro-movements of one edge 
with respect to the other during function, such as walking, to be successful. The stability of 
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device depends on its shape, and material, and on the number, position and inclination of the pins 
that secure the device to the tibia. The best position and inclination of pins for the stability of the 
device depends on the spot chosen, that is on the local property of the tibia. The anisotropy and 
non-homogeneity of the tibia make a difference with respect to the screw loosening while 
walking. In fact, the screw may or may not get loose if the chosen spot is more or less resistant to 
the force that it takes for the pin to get into place, if one inclination is chosen instead of another 
one, if the spot is prestressed in one direction or another. The question of osteotomy fixation 
stability cannot be fully studied with a computer bone model that does not take into account the 
hierarchical structure of bone that renders bone anisotropic and non-homogeneous. Because if 
the mode] assumes less than that, the local information is lost and the bone shows the same 
properties where it should not. 

Another example refers to cemented implants. The local bone conditions affect the bone- 
implant interface. The loosening rates in cemented implants, especially in younger, active 
persons, is partially due to the local bone mechanical properties. This problem has led many 
investigators to pursue methods of cementless fixation. In the meantime a great deal of attention 
is being focused on the bone-implant interface and the factors affecting its strength. A thorough 
solution to the problems involves the knowledge of the local bone mechanics. 

Simulation of Fracture Propagation 

The fracture propagation model of either compact or cancellous bone under specific 
loading follows the same steps as the fracture propagation model of single osteon samples under 
torsion (See, Example 2, Part E, Steps 1-19). The computer program may be based on any 
suitable simulation program, for example, a Monte Carlo simulation. The fracture propagation 
steps are applied to the finite element mesh for the compact or cancellous bone in question, 
instead of to the finite element mesh of single osteon samples. 

The purpose of the fracture model isto show that cumulative micro-cracking, de-bonding, 
void growth and fiber breakage associated with repeated loading of osteons causes a progressive 
loss of stiffness and pinching, and increase of energy absorption. 

The fracture model reflects the following hystological/physiological observations. Fluids 
occupy vascular canals, canaliculae and lacunae, which are interconnected. The flow of liquids 
under stress can absorb large amounts of energy, increasing the toughness of bone. Large strains 
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may be accommodated by the organic phase (e.g. collagen, mucopolysaccarides). When a strain 
is sufficient to cause cracking, the organic phase may also contribute to the dissipation ofenergy 
at the front of a propagating crack. Crack propagation also appears to be arrested in the presence 
of canaliculae and lacunae. In fact, when the crack gets to a hollow space, it just stops because at 
the hollow space there is no more resistance, no more material to rip. Therefore, discontinuities 
to some extent increase the robustness of bone rather than increase its tendency for brittle fracture 
(Currey, 1 962). In the case where a crack enters a discontinuity, the front tends to be blunted, 
hence reducing the stress concentration factor (i.e.- the level of stress necessary to create a crack) 
and slowing crack propagation. When a crack is forced to enter a vascular canal, the radius at the 
tip of the crack becomes larger. Lacunae are probably more likely to act as stress concentrators 
than canaliculae because of their generally ellipsoidal cross-section and because they are 
generally oriented normal to the long axis. Stress concentrators are define as entities that raise 
the stress concentration factor. However, their much smaller size precludes them from acting as 
fracture initiators (i.e. causes for the structure to begin fracturing) until or unless plastic 
deformation has created cracks at the tip. Fractures spread along cement lines and lamellar 
interfaces. 

Additional Applications 

The model also can be used with complimentary applications and technologies. An 
example includes, but is not limited, to combination with software to model soft tissue (such as 
the one developed by the company lnfocus, Sylicon Valley, California). Material analogs of bone 
can be obtained by means of 3D printers (see e.g., the printers manufactured by Stratasys- 
SDprinting in Eden Prairie, Minnesota). Implants and distraction devices will be manufactured 
by computer guided robots. See, e.g., Mah and Hatcher, 2000. The present model will provide 
the model of the bone structure to be distracted and during distraction. 

The application of modeling to imaging (e.g., clinical MRI and CAT scan x-ray imaging) 
of human bone offers the prospect of a qualitative leap in the predictability, effectiveness, and 
convenience of surgical, orthodontic, orthopedic and other medical interventions. Embodiments 
of the mode! can enable medical professionals, based on patient specific data, to visualize how 
bone in various parts of the body will grow and heal in response to medical interventions. First, 
the current lack ofbone local mechanical properties impairs the comparison between natural and 
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synthetic bones. Second, the current lack of knowledge of mechanical properties, strain and 
stress distributions throughout the bone, impairs the research for new synthetic bones to move 
towards the same properties. For instance, the latest synthetic long bone is made out of fiber- 
reinforced glass (Szivek, 2000) of unknown local mechanical properties. There are no reports of 
synthetic porous structures with the interconnecting pores having the same stiffness and strength 
characteristic of human trabecular bone. Even the most popular synthetic closed-cell 
polyurethane foams (such as Daro, Butler, Wl) which have a structure that shows similarities 
with human trabecular bone are homogeneous in theory and with inhomogeneities difficult to 
control in practice. In any event the non-homogeneous hierarchical structure of human bone is 
not even close to being imitated. Third, the current lack of knowledge of mechanical properties 
and of strain and stress distributions throughout the bone, impairs bone reconstruction, bone 
grafting, placement of screws, insertion of prostheses. 

The invention is applicable to the bones of other vertebrates whose bone structure 
somewhat differs from that of humans. For instance, the invention would give valuable results 
on the prevention and healing of fracture in equine bone. Currently, the micromechanical bone 
studies of vertebrates are scarce, often limited to a few small animals, such as mouse, dog, and 
sheep. Because of that, the results on human bone microstructure are erroneously used in studies 
of vertebrates to which they do not apply (Riggs et al, 1993a, 1993b). 

The present invention provides a more realistic prediction of the macroscopic bone 
mechanical properties, strain, and stress distribution than computer models based on omission of 
either anisotropy or non-homogeneity of bone. Moreover, this invention provides more realistic 
prediction than purely mathematical models, that is models based on hypotheses, which are not 
based on experimentation. The literature is full of research on bone microstructure, which 
employs purely mathematical models of osteon behavior (Pidaparti and Burr, 1992). Such 
approach is limited, often unrealistic and does not always predict biological phenomena. The 
invention is flexible so as to include new experimental findings of bone structural and 
mechanical properties. This ensures the invention's realistic characteristic insertion of 
prostheses, etc. 

EXAMPLES 
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The present invention will be better understood by reference to the following 
Examples, which are provided by way of exemplification and not by way of limitation. 
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EXAMPLE 1 

To produce a model of the present invention, compact bone is subjected to any method that 
may produce non-invasive slices of biological structures that are known in the art (/>., jiCT-scan or 
micro computerized tomography). Images then are stored in a computer and a 3-dimensionaI- 
reconstruction is applied using a standard method known in the art (see e.g. Materialise, XYZ 
Scientific Applications, Inc. Livermore, California). The high resolution of a |aCT-scan (about 30 
jam) allows for determination of the outline of osteons, of osteons 3 vascular canals and interstitial 
bone. Also the 3D-reconstruetion shows varying shades of gray, which represent the degree of 
calcification. Osteons are filled with structure by means of the two lamellar types (bright and extinct 
lamellae), which have been previously assembled. The criteria by which the lamellar structure is 
drawn into each osteon follows the distribution of alternate osteons (Hert, et al 1994) and the 
distribution of dominant collagen fibril directions (Portigliatti-Barbos, 1983, 1984, and 1987; Boyde 
elaU 1984; Ascenzi A. ei aL, 1987a and 1987b; Ascenzi A., 1988; Carando etaL, 1989 and 1991). 

The structure of the osteonic lamellar model consists of a laminate whose length, width and 
height correspond to a cylindrical shell circumference, thickness, and height (Fig. 9a). The structure 
of the lamella within the interstitial bone is modeled as a portion of the osteonic lamellar model (Fig. 
9b). The layers are unidirectional fiber-reinforced laminae (Fig. 9c) of the same matrix and fibers. 
The matrix and fibers, Le. individual components of the hierarchical structure, and not the 
microstructure as a whole, are each treated as homogeneous and isotropic. 

All fibers, of which there are two types, are assumed to be circular in cross-section, randomly 
distributed in the transverse plane. The fibers of the first type have a diameter of about 800 A and 
fibers of the second type have a diameter of about 200 A. The fibers of the first type are proposed to 
be perfectly embedded in matrix (they are idealized to have essentially no gaps between them and do 
not move relative to each other). The fibers of the second type are proposed to be perfectly 
embedded in matrix only when bone undergoes physiological static loading. When bone undergoes 
physiological dynamic loading the fibers of the second type are given the option to move with 
respect to matrix in which they are embedded. Such displacement is specified following 
experimentation, such as boundary condition or de-bonding experiments. The experimentation may 
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also dictate further conditions for the relative position of the two fiber types. Examples of such 
experiments are discussed below. 

In the present model, the lamina with fiber inclination y is named y-lamina. The thickness of 
the dark lamella ranges between 7 and 12 um (Ascenzi A. el al, 2000). It is described by the 
sequence [82, -82] (Frasca et al , 1 977). The notation [82, -82] refers to two y-laminae where y=82, - 
82. The thickness of the bright lamella ranges between 4 and 7 urn (Ascenzi A., et al, 2000). It is 
described by the sequence [-61 .5, -41,-20.5, 0, 20.5, 41, 61 .5] (Ascenzi M.-G., 1999b). 

For the matrix, Young's modulus of 1 1 4 GPa, Poisson's ratio of 0.27, and ultimate strength of 
0.59 GPa are assumed for hydroxyapatite (Katz and Ukraincik, 1971). For the fibers of the first type, 
Young's modulus of 1.2 GPa, Poisson's ratio of 0.35, yield strength of 0.002 GPa are used for 
- collagen (Currey, 1 969). For the fibers of the second type, Young's modulus of 1 . 1 GPa, Poisson's 
ratio of 0.23, are used for mucopolisaccarydes (Bourne, 1971). Depending on the degree of 
calcification, the matrix occupies up to 40% of the lamina volume without voids (Bonfield and Li, 
1 967). The cement line is modeled as homogeneous and isotropic: the Young's modulus of 70 GPa 
and Poisson's ratio of 0.27 (Philipson, 1965; Schaffler et al, 1987). 

Since osteons vary with respect to the distribution of dark and bright lamellae, the model of 
an osteon with a specific distribution of dark and bright lamellae is obtained by assembling the 
model of dark and bright lamellae so as to follow the osteon's particular distribution. For instance, a 
model of the longitudinal osteon, which consists of dark lamellae, is made of 12 laminae. The fiber 
inclination angle changes from 82 ° and - 82° six consecutive times. A model of the alternate osteon, 
which consists of alternating dark and bright lamellae, is made of 36 laminae. The fiber inclination 
angle increases by 20.5 0 from -82 0 to 82 0 and then decreases by 20.5 0 from 82 0 and - 82° four 
consecutive times. 

Information included for the present model may not be currently available for all bones that 
are evaluated. Any information that is needed for the practice of this invention may be obtained by 
experimentation using methods that are standard in the art. Additionally, methods that may be used 
to evaluate bone in one species may be used to evaluate a similar bone structure in another species. 
For example, the distribution of dominant collagen bundle directions is available for the shaftof the 
human long bone but not for other vertebrates nor for the mandible. For any compact bone the 
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distribution of dominant collagen fibril directions can be obtained by applying the method of Boyde 
etui (1984). 

For cancellous bone the same method can be applied after embedding (soaking and letting 
dry) the bone in a conventional resin used for the specimens examined under the electron 
5 microscope. Such resin should not change the microscopic characteristics (birefringence) of the 
specimens, so that the image of the collagen bundle and hydroxyapatite needle directions under the 
polarizing microscope is not altered by artifiacts. An example of such a resin includes, but is not 
limited to, epoxy. Note that application of the invention to cancellous bone will model lamellae that 
form trabeculae, as compared with the osteons of compact bone, however trabeculae and osteons can 
10 both be modeled in terms of lamellae. 

From the above-mentioned mechanical properties of matrix and fibers {e.g. t Young's modulus 
and PoissoiVs ratio) the same types of mechanical properties for lamellae under various load types 
(such as tension, compression, shear, and torsion) will be deduced by means of standard fiber- 
reinforced laminate methods known in the art (see e.g., Jones, 1975; Vinson, 1993; Antman, 1995). 

o 

15 Based on mechanical properties of the lamellae, homogenization theory will allow for the 

deduction of the osteon, osteon group, and interstitial bone mechanical properties for compact bone 
and trabecular mechanical properties for cancellous bone. The mathematically computed mechanical 
properties of lamellae, osteons, osteon groups, interstitial bone, and trabeculae are compared to the 
experimental results. If the experimental results are not available for the particular bone to which the 

20 invention needs to be applied, the properties may be determined using the methods for mechanical 
testing as described herein. The mechanical properties of lamellae, osteons, osteon groups, 
interstitial bone, and trabeculae are used as input for the homogenization methods to deduce the 
mechanical properties of the desired macroscopic bone. 

Results are included of a finite element model, which allows for the assessment of the 

25 mechanical properties of the sample. The sample dimensions before and after testing allow for the 
formulation of an equation that describes the deformation from the shape before testing to the shape 
after testing. The deformation equation allows for the computation of the strain distribution 
throughout the sample. For example equations see Antman, 1995. The combination of such strain 
distribution with the experimental diagrams, the known sample structure before testing, and the 
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fracture patterns after testing allow for the computation of the elastic properties through standard 
finite element methods. Statistical student t-lest (Moore and McCabe ; 1989) is run across the 
sample's results to allow for comparison of mechanical properties across the samples and to allow 
statistical conclusions. 

5 These studies provide the mechanical properties of all the hierarchical orders. Therefore, the 

mechanical property distribution throughout the bone in terms of the microstructural components is 
known. The finite element method is applied (see e.g. the software package Abaqus) to compute the 
bone response to any given force acting on it. Boundary conditions are entered as assumptions into 
the finite element method. The first step is to create a 3-dimensional mesh (see e.g. Couteau etai, 
10 2000). 

The bone overall shape is filled with "elements". These elements are used to represent the 
osteons present in the bone. For example, a hollow cylindrical portion of an osteon with an inner 
and diameter of 40 u.m, an outer diameter of 220 urn, and height of 500 um, is filled with about 
600,000 elements. Mechanical properties and boundary conditions are the method's input. 

1 5 Boundary conditions express the movements of the boundaries of the various hierarchical structures 
under loading. For example, dynamic loading evidences bone's viscous behavior. The literature 
points to mucopolysaccarides or perhaps collagen as the microstructural component responsible for 
the viscosity. In the structural part of the invention the second type of fibers models the 
mucopolysaccarides. The fibers of the second type are free to move. Such movement at the 

>0 interface between the fibers of the second type and the matrix is expressed by a boundary condition 
(to be determined experimentally). Another example, the behavior of the cement line under loading 
r is the boundary conditions for the osteon and the interstitial bone between which it lies. If the 
boundary conditions of a specific bone, to which one wants to apply the invention, are not available 
in the literature; they can be assessed experimentally by applying methods described herein or that 

5 are well known in the art. The software application gives as output the strain and stress distributions 
throughout the bone. 

The mechanical properties of compact bone microstructure (lamellae, single osteons, osteon 
groups, single trabeculae) can all be experimentally found with the following method (other 
acceptable methods, including non-invasive methods are available in the literature). Human 
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cadaveric bone aged between 20-50 is obtained according to the US regulations. The cadaveric 
bones are chosen either free of pathology to apply the invention to normal bone or with a specific 
pathology to apply the invention to a specific pathology. The bone marrow is removed by standard 
anatomical techniques (Wickramasinghe, 1975). At least 15 samples of any of such structures 
5 (lamellae, single osteons, osteon groups, single trabeculae) are isolated from the surrounding bone. 
The samples have about the same size and shape. The shape is a parallelepiped, a cylinder, or a 
hollow cylinder (depending on the chosen structure) with lugs (see e.g. Fig. 4c) for mechanical 
testing. Sample preparation and selection of compact bone microstructure is achieved by the 
methods of Ascenzi A. et al (1 994, 2000). For example, although any technique can be used, the 
1 0 method of Ascenzi A. et al. (1 994) is preferably used to isolate osteons. The preferred form chosen 
for isolation of osteon samples is a cylindrical shape around the vascular canal. In general, the shape 
and location of a structural sample are chosen in such a way so that all the properties of the structure 
are preserved. Mechanical testing of osteon samples (Ascenzi M.-G. et al, 2000) may include, but is 
not limited to, monotonic and cyclic testing in tension, compression, shear, bending and torsion. The 
1 5 methods conducted as described in Ascenzi A. and Bonucci, 1 967; Ascenzi A. and Bonucci, 1 968, 
1972; Ascenzi A. et al, 1990, 1994; Ascenzi A. et al, 1985, 1997, 1998 have proven themselves 
successful. The testing is conducted within the elastic range and beyond the elastic range to study 
fractures. Sample preparation and selection of cancellous bone (single trabecula and trabecular 
groups) is achieved by any of the methods whose bibliography is listed in Mente, 2000. Each sample 
20 is measured (the three dimensions for the parallelepiped; base radius and height for the cylinder base 
inner and outer radii and height for the hollow cylinder) before and after isolation and before 
mechanical test. Change in dimensions before and after isolation and before mechanical test shows 
existence of prestress. The structure of the sample is assessed before or after the mechanical testing 
(Ascenzi M.-G. et al , 2000). The samples are tested mechanically under physiological conditions, 
25 that is wet at 21° C. Since both compact and cancellous bone are viscoelastic, the results of 
mechanical testing are time-dependent (Sasaki, 2000). Consequently the strain rate and testing 
frequency need to be prechosen and the computer modeling depends on such choices. The stress- 
strain experimental curves (either monotonic or cyclic) through the elastic arid plastic ranges are 
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evaluated and recorded. After the mechanical test, the bone samples are measured and observed 
under the optical microscope for fracture patterns. 

A trapezoid is cut around each osteon sample (see Figure 13). For osteon immobilization 
during lamellar isolation, a portion of the bone material inside the trapezoid around the osteon is 
glued with Kemi®Cyak adhesive to a slide. The bright and the extinct lamella at the periphery of 
each osteon are dissected with a razor-sharp microscopic blade, obtained by filing a steel needle. To 
avoid fracture formation during straightening of each lamellar sample, such operation is performed 
gently on wet samples while checking under an optical microscope. The selection of external 
lamellae, of lesser curvature than internal lamellae, decreases the risk of fracture formation during 
flattening. The ends of each flattened sample are secured to two supports. The samples are 
* measured as previously described and examined under an optical microscope to assess defects. 

The mechanical loading on lamella on wet samples is conducted statically at 21°C to 
complete rupture, with current model of the microwave extensimeter. Flattened bright lamellar 
sample is expected to resist tension along its length better than extinct lamellar sample. Indeed, the 
bright lamella is hypothesized to contain collagen bundles transverse to the longitudinal lamellar axis 
when enclosed in bone. Therefore, the enclosed lamella's transverse bundles strengthen the flattened 
lamellar sample in the direction of its length. The extinct lamella is hypothesized to contain collagen 
bundles parallel to the longitudinal lamellar axis when enclosed in bone. Hence, the enclosed 
lamella's longitudinal bundles after flattening are a source of sample weakness in the direction of its 
length because they are transverse to the loading direction. Fracture patterns of ruptured samples are 
studied under an optical microscope. Observation of fractures in ruptured samples will allow 
- formulation of hypotheses on fracture nucleation and growth. 

Since the interest is to test the isolated and flatten lamellar samples mechanically, the stresses 
present in the flatten lamellar samples before testing are assessed. To compute the stresses in the wet 
flattened lamellar samples, a computerized geometric-material model of a bright wet lamellar sample 
and of an extinct wet lamellar sample will be constructed, separately before and after isolation and 
flattening. The bright lamella includes prestress. It is hypothesized that the the stresses in the flat 
bright lamella are larger than the ones in the extinct lamella. Additionally, by taking into account 
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that the periosteous is prestressed in tension, it may very well be that the outer circumferential 
system is prestressed in tension, too. 

The geometry of the model is based on (1) dimensions (inner and outer radii, height, and 
dimension variations) of wet lamellar samples before isolation from surrounding alternate osteon and 
5 after isolation and flattening (width and length., and dimension variations) and (2) structure of 
lamellar sample. Therefore, dimensional measurements are needed. The structure model also is 
based on the lamellar structural components' arrangement. Therefore, lamellar structure under a 
confocal microscope will be assessed. 

10 EXAMPLE 2 

A. Sample Preparation, Measurements, and Experimentation 

Sample preparation and selection in this example will be achieved by applying the Ascenzi 
k. etal. (1994) methodology. Accurate sample preparation is important, and the Ascenzi method is 
1 5 preferred. The femoral shafts of human corpses free from evident skeletal faults will provide the 
bone material for the study. Longitudinal shaft segments about 30 mm long will be first sawn off, 
and longitudinal sections slightly thicker than an osteon (350 urn) will then be prepared using aLeitz 
rotating-saw microtome. A continuous water-spout will be incorporated to prevent any overheating 
of the material. Osteon samples will then be isolated from the sections. The features determining 
20 sample selection will be the degree of calcification and the orientation of the collagen bundles and 
crystallites. Microradiographic examination, preferably according to Amprino and Engstrom (1 952), 
will allow the selection of fully calcified osteon samples. Two types of osteons will be selected. 
They correspond to two different collagen bundle patterns in fiber orientation in successive lamellae. 
Under the polarizing microscope, one type, the longitudinal osteon, has a predominantly dark 
appearance in cross section (Fig. 4a); whereas the other type, the alternate osteon, reveals alternately 
bright and dark lamellae (Fig. 4b). 

When bone sections are cut longitudinally, the two osteon types are easy to recognize 
provided that the thickness of the sections is much less than the diameter of an osteon. Longitudinal 
osteons appear to be almost uniformly bright under the polarizing microscope, while alternate 
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osteons show alternatively bright and dark lamellae. When the thickness of the bone section differs 
little from the mean diameter of the osteon, concentric lamellae overlap, thereby reducing or 
precluding the visibility of dark lamellae, and leaving open the possibility that an alternate osteon 
may have a bright appearance. As a result identification only becomes certain once a cross section 
has been cut from the osteon using a microscopic drill (Ascenzi A. et al., 1994). Hence, cyclic 
loading must be performed before undertaking positive identification of the osteon type. 

With reference to the position and orientation of the haversian canal, it is necessary that the 
canal lie midway between the surface of the cylindrical sample and parallel to it, so that torsion is 
applied around the osteon axis. This calls for the preliminary separation of the osteon sample, e.g. 
by application of a technique described below. This technique allows the position and orientation of 
the canal to be calculated by measuring its distance from the outer surface of the sample at various 
levels and rotational angles. 

The samples are isolated in two stages. During the first stage the sample, consisting of the 
central portion of an osteon, 500 urn in length, with the ends penetrating into two rectangular lugs, is 
separated from the bone section using a device as described in Ascenzi A. and Bonucci (1968) and 
Ascenzi A. (1990). As isolation of the central portion of the osteon is achieved by drilling, its 
section has a coarse, square shape. During the second stage, a micro-grinding lathe is used to give 
the central portion a cylindrical form, with the haversian canal running through it axially. The lathe 
to be used was designed and developed by the CECOM Company and is described by Ascenzi A. et 
al. (1994). The device grinds the sample by a minute steel blade whose edge, 500 um long, is equal 
to the length of a coarsely isolated sample. The forward and backward movements of the blade are 
monitored by a micrometer. The length and other dimensions of the various samples were kept 
virtually constant; one criterion for the choice of the samples is that their haversian canal measures 
40+3 um in diameter. Additionally, a stopper controls the forward and backward movement of the 
steel blade on the micro-grinding lathe to provide a series of samples whose external diameter equals 
21 0±3 u.m. This provides a precise comparison of samples' torsional properties (Ascenzi A. et al, 
1990). Osteons are not uniform in dimensions. With the dimensions carefully controlled and 
standardized to exclude defects and other structures, the material rather than structural properties are 
determined for the osteons. This information can then be applied to osteonal structures of varying 
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dimensions under the assumption of homogeneity at the level of the osteon rather than for the 
macroscopic specimen. 

The relative dimensions of the osteon samples may not appear to conform to those 
conventionally suggested for material testing. They reflect conditions made necessary by the 
distinctive nature of bone microstructure. In particular, 500 urn is the maximum length compatible 
with the avoidance of Volkmann's canals in the wall of the specimen. An external diameter of 210 
Urn is the maximum dimension possible that ensures that portions of the neighboring structures are 
not included in the sample as a result of irregularities in the thickness of an osteon. The internal 
diameter of fully calcified osteons averages 40 um. 

Fig. 4c shows a completely isolated osteon sample held within rectangular lugs. The lugs 
allow the sample to be firmly attached to the device while hysteretic loops are recorded. The central 
portion of each sample will be only 500 um long; consequently, the sample will not include 
Volkmann's canals which would behave as discontinuities. In addition, the osteon sample selection 
criteria includes that the vascular canal should run strictly parallel to and equidistant from the surface 
of the cylindrical sample and that there should not be small surface defects. The canal's position and 
orientation are assessed by checking the distance between vascular canal and external surface of 
sample at various rotational angles and levels. To exclude the presence of small surface defects that 
could alter the shear modulus values in torsional testing, each sample is subjected to careful optical 
microscope examination. Severe criteria are set for osteon sample selection. Osteon types can only 
be identified from a prepared cross section only after a sample has been tested. This means that to 
obtain 60 samples divided between those containing longitudinal and alternate osteons, which 
satisfactorily complete the procedures adopted for the recording of the hysteretic loops under torsion, 
it will be necessary to prepare between 800 and 1,000 samples. 

The apparatus is an adaptation of the device for testing osteons under torsion to failure 
described in Ascenzi A." el a!., 1994), and further described in Figures 10 and 11. This device 
consists of a rotational axis, point (1) in Fig. 1 1, with two sets of jaws, point (2) in Fig. 1 1, which 
grip the specimen during testing. The jaws are oriented along the same axis but none of them are 
free to move axially. This sets up an axial loading effect, which could influence absolute 
measurements but may be neglected when, as in this investigation, comparative measurements are 
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considered. 

One set of jaws is fixed, while the other turns in synchrony with a wheel, point (4) in Fig. 1 ! , 
measuring 61 mm in diameter. In order to minimize the rotating friction of the turning jaw, a 
pendulum loading system is adopted. The axis of the pendulum loadmg system is indicated as point 
(5) in Fig. 1 1 . The frictionless fulcrum of the pendulum loading system is the tip of a hard metal 
wedge, point (3) in Fig. 1 1 . The maximum oscillation of the pendulum is fixed at 55° A tungsten 
thread, whose section measures 20 urn in diameter, winds around the rim of the wheel as a series of 
0. 1 gram weights is attached i ncrementally at one end of the tungsten thread. Weights are added one 
by one until failure occurs for monotonic torsional loading. The load limit is chosen so that the 
corresponding torque is approximately equal to the middle value between the maximum possible 
elastic torque and fully plastic torque. Preliminary trials indicate a load limit of 0.9-1.0 gram for 
fully calcified longitudinal osteons and of 0.8-0.9 gram from fully calcified alternate osteons and of 
0.6-0.7 gram for decalcified longitudinal osteons and of 0.5-0.6 gram from decalcified alternate 
osteons. Once the load limit is reached, the weights will be detached one by one. The procedure will 
then be repeated at the other end of the tungsten thread. In this way, the osteon specimen rigidly 
clamped at one end is progressively twisted at the other end by a torque in counterclockwise and 
clockwise directions alternatively, so as to achieving cyclic loading. The interval between the 
application of two consecutive weights will be kept constant at 4 sec. A stereoscopic microscope 
will be used to verify that the axis around which torsion occurs coincides with the osteon axis. The 
aim of this operation is to check that the center of each jaw corresponds to one end of the osteon 
sample canal. The angle through which one end of the specimen twists relative to the other during 
testing is measured by applying an optical method based on the reflection of a laser beam from a 
small mirror attached to the rotating set of jaws. The variations in the angle of twist are read on a 
graduated scale placed 1 60 cm from the device. The precision and accuracy of the graduated scale 
coincide with those of the apparatus, as checked by applying experimental procedures. Because the 
diagrams obtained when testing begins in the counterclockwise (or positive) direction should look 
essentially like the diagrams obtained when testing begins in the clockwise (or negative) direction, 
all the diagrams will be recorded starting in the counterclockwise direction, according to the standard 
practice reported in the literature. 
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The cycles applied to each specimen will vary in number and they will be interrupted before 
the spontaneous sample rupture. Preliminary specimen testing taken to rupture gives an indication of 
the angle-of-twist values at rupture and therefore give an indication to stop cycling when angle-of- 
twist values get close to the preliminary rupture values. The final cycle will therefore be interpreted 
as the cycle preceding the interruption of the experiment; in consequence, it has no physical 
meaning. Issues related to osteon fatigue-life will not be part of this study because they require the 
ultimate destruction of the specimen which would prevent the ability to properly identify the sample 
after cyclic loading. Interruption of the experiment prior to rupture is necessary because, as 
previously noted, the osteon type can only be identified with certainty by preparation of a cross- 
section before rupture, but after testing. Before preparation of the cross-section for identification, 
samples will be examined under the optical microscope to analyze the nature and size of lesions. 
Osteon samples are examined under an optical microscope as whole, even though to do so samples 
need to be removed from the torsion device; such removal can cause structural changes, e.g. partial 
closure of cracks. The nature of fatigue-damage cannot be verified directly under microscope. In 
fact, fractures due to cycling and those due to osteon sawing to produce the section are not 
distinguishable. There is no universally accepted technique available to allow qualitative 
observation of fatigue-damage. The only experimental alternative is x-ray diffraction, which 
provides only quantitative indication of fatigue-damage (Ascenzi A. etal, 1998). 

Of 60 osteon samples, 7 longitudinal and 7 alternate osteon samples will be decalcified by 
treatment in a versene solution buffered to pH 7 (Ascenzi A. et al. , 1 967). Measuring the increase in 
birefringence at regular intervals of time will check the decalcification (Ascenzi A. and Bonucci, 
1964). After hydration of the material with saline solution, 20 longitudinal and 20 alternate fully 
calcified osteon samples, and the 7 longitudinal and 7 alternate decalcified osteon samples will be 
tested wet as described above. The remaining 6 osteons will be tested on the first and second cycles 
only and their dimensions will be measured (by means of a micrometer under an optical microscope) 
after cycling and before the sectional cutting for osteon type determination. The mechanical testing 
will be conducted at 20°C. The specimens will be maintained wet during testing by continuous use 
of a micropipette. 

Upon completion of the experimental portion of the research, hysteretic values of torque vs. 



WO 02/33679 



PCT/US01/32474 



angle of twist will be plotted for each osteon sample. 

The proposed mechanical testing is performed quasi-siatically so rate dependencies are not 
expected in the material response. To confirm this, a series of twist and hold experiments will be 
performed as a preliminary test to see if significant strain relaxation (creep) occurs in the time frame 
that it takes to complete an experimental cycle. The presence of creep would manifest itself by 
maintenance of the shape of the curve but with a clear translation of the entire hysteresis loop. If 
creep is not present, the area under the hysteresis loop is the energy absorption. If creep is present, it 
will be accounted for by the addition of viscoelastic terms. 

B. Mathematical Analysis of Hysteretic Curves 

To establish analogies and differences among plots, the characteristics of torque vs. angle of 
twist plots are quantified by means of a polynomial approximation for each half-cycle. (Shiga et al., 
1970; Ascenzi A. etal, 1997). 

Measurements and Plots 

(a) Least-square regression will be applied to the data for each half-cycle to identify the best 
polynomial approximation of at least second order. The degree cannot be 1 because the half-cycle 
involves the inelastic range which is characteristically nonlinear. 

(b) The goodness of the approximation will be determined through analysis of the residuals 
and computation of the percent variation in torque explained by the regression. A lower bound for r 2 
is set to 0.98. The degree of the polynomial approximation at (a) will be increased in a stepwise 
fashion until a good-approximation is found for all half-cycles. 

(c) Let n be the smallest integer for which there is a good-approximating polynomial for 

n 

all half-cycles. Let the polynomial equation be y(x)= Z a r x f . Note that such an equation does 

i=o 

not include symbolism that denotes the half-cycle; this is done to add clarity and does not lead to 
confusion. A literature search suggests that n might equal 3 or 4 (Shiga et al, 1970 and Ascenzi et 
al, 1997a). Each coefficient a, for any given cycle will be plotted with respect to the maximum 
angle-of-hvist on that cycle, for negative and positive torque half-cycles separately for each osteon 
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type, to visualize the value's variation as the number of cycles increases. 

(d) The first derivative of the polynomial, y'(x)- S ia, x (i-, \ will be taken to represent the 

i=o 

stiffness of each individual half cycle. 

(e) Compute relative extrema of y'(x) on each half-cycle. Because of increasing structural 
damage, stiffness is generally either always decreasing on both negative and positive torque half- 
cycles or it shows a sharp change from decreasing to increasing only once on each half-cycle. 
Therefore, y l (x) should show no relative maximum and a marked relative minimum. If such relative 
minimum exists, the graph of y(x) shows an inflection point, i.e. pinching is present. Let h m and k m 
denote twist and torque values at the stiffness relative minimum of the half-cycle. In these notations, 
y'(h m ) denotes the minimum value of stiffness. 

(f) Compute y'(x b ) and y(x c ), the stiffness at the begin and at the end of each cycle. 

(g) Compute the absolute maximum of the function |y rr (x)|/(R(y , (x))2)3/2 on any given half- 
cycle. This is the maximum curvature per half-cycle and will be denoted by m c . 

(h) Compute the difference between the integral of y(x) over the first half-cycle and the 
integral of y(x) over the second half-cycle. Such value approximates the area of the region bounded 
by the first cycle. This is the energy absorption during the first cycle and will be denoted e,. 

Differences in hysteretic behavior between longitudinal and alternate osteons will be 
established by studying the distributions of maximum twist, h m , k m , y'(h m ), y\x b yyXx e % rn C) and e t 
computed above. The statistical t-test, paired or unpaired, will be applied on means of distribution or 
on the mean of the distribution's logarithm if the distribution is not normal. 

The mean of the differences of twist limits between the last and first cycle obtained from the 
experimental diagram will be computed. It will be compared between negative and positive torque 
half-cycles for longitudinal and for alternate osteons, separately. The magnitude of such value should 
be smaller for longitudinal osteons than alternating osteons because longitudinal osteons resist 
torsion better than alternate osteons. 

The signs of h m (and k m? respectively) at the first and last cycle will be analyzed. h m (and k m , 
respectively) should have the same sign for the two osteon type separately, up to possibly a few 
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samples. This would indicate that the twist at minimum stiffness should not change much at all 
within ail negative and positive torque half-cycles separately, for both osteon types. This would be 
in agreement with the tight and well organized osteon structure. 

A paired two-sample t-test will be applied to the means or the values of h m (and k„„ 
respectively) at each of the cycles in any given set of corresponding cycles of the two osteon types. 
This compares the values of minimum twist (and torque, respectively) for negative and positive 
torque half-cycles. 

A paired two-sample t-test will be applied to the means of y'(h m ) and the coefficients of y' (x) 
for the two osteon types, separately at each of the cycles in any given set of corresponding cycles of 
hysteretic diagrams. This compares the minimum value of the stiffness between the negative and 
positive torque half-cycles. 

An unpaired two-sample t-test will be applied to the means of h m (and k m , respectively) for 
the two osteon types at each of the cycles in any given set of corresponding cycles of hysteretic 
diagrams. This will compare the values of twist (and torque, respectively) at the inflection point for 
the two osteon types. 

A unpaired two-sample t-test will be applied to the means of y'(h b )-y'(h e ) for the two osteon 
types at each of the cycles. This will compare the stiffness decrease within a given cycle between the 
two osteon types. 

A paired t-test will be used on the mean of the difference of m c for negative and positive 
torque half-cycles at each of the cycles in any given set of corresponding cycles of hysteretic 
diagrams. This will compare the maximum value of the curvature of the stress-twist diagram 
between negative and positive torque half-cycles, for the two osteon types, separately. The 
maximum value of curvature of the torque-twist diagram is expected to be larger on the positive 
than negative torque half-cycle for longitudinal osteons; whereas, the maximum values of the 
curvature of the torque-twist diagram under negative and positive torque half-cycle for alternate 
osteons should show no difference. As a consequence of the reduced stiffness the energy absorption 
should be larger for longitudinal than for alternate osteons. 

A paired t-test will be used on the difference of the means of each coefficient a; of y' at last 
and First cycles, on both negative and positive torque half-cycles, separately. For both osteon types, 
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the value of y' is expected to decrease from the first to the last cycle on both negative and positive 
torque half-cycles for any value of the twist x. This test will measure stiffness degradation here 
defined as the decreasing of stiffness at any given twist value on either a negative or positive torque 
half-cycle as the number of cycles increases. 

The existence of the value h m shows the S-shape of the half-cycles that identifies pinching. 
Pinching is expected to be present for each cycle for both types of osteon. If pinching is present, 
pinching degradation will be computed. 

A paired two-sample t-test will be applied to the means of y'(h ra ) for the two half-cycles of 
any given cycle. This compares pinching degradation for the two osteon types separately. Pinching 
degradation at any cycle is the reduction in stiffness from its value at the deflection point of the 
negative torque half-cycle to a lesser value at the deflection point of the positive torque half-cycle. 

An unpaired two-sample t-test on the means of the value of minimum stiffness y'(h,„) for 
longitudinal and alternate osteons at any given half cycle. This compares pinching degradation 
between the two osteon types. 

An unpaired two-sample t-test will be applied to the means of e, for longitudinal and 
alternate osteons. This compares energy absorption between the two osteon types at each of the 
cycles in any given set of corresponding cycles of hysteretic diagrams. 

The mechanical meaning of some of the parameters used in the above analysis of 
experimental diagrams (e.g. stiffness, energy absorption) will be made clear as such parameters will 
be correlated to ultrastructural behavior during fracture propagation in the model described at section 
IV. Up to this point they are comparative measures of behavior between longitudinal and alternate 
osteons under torsional loading. 

C. Interpretation 

A structural and biological interpretation of the shape of torsional hysteretic loops in osteons 
through the results of the previous steps uses a segment representation of each cycle of the 
curvilinear recorded diagram. 

If, as it is anticipated, pinching exists, the bilinear model of Fig. 12 is appropriate. Here 
points B and E approximate endpoints of the negative torque half-cycle, while points E and H 
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approximate endpoints of the positive torque half-cycle; segments DC and FG approximate tangent 
lines to the curves at the inflection points. The three segments modeling the twist decreasing branch 
show that stiffness decreases (along segments BC and CU) to a minimum value and then increases 
(along segments UD and DE). Similarly, the three segments modeling the positive torque half-cycle 
show that stiffness decreases (along segments EF and FL) to a minimum value and then increases 
. (along segments LG and GH). 

The slope of segment DC on the negative (FG on the positive, respectively) torque side of the 
bilinear model is smaller than both the slopes of segments ED and CB (EF and GH, respectively) and 
therefore responsible for a contraction of the cycle that constitutes pinching. The existence of 
pinching resides in the torque-angle-of-twist branch AB of the primary curve, where lesions appear 
, _ as a result of yielding of components of the bone matrix under load as the angle-of twist increases to 
the right. Reversal of loading is required to close the lesions; this will occur once the minimum 
■angle-of-twist of the unloading portion BC of the curve is exceeded. Once point C is passed, 
stiffness, as negative torque is decreasing, shows a progressive, slight, unsteady increase to point D. 
The lesions are then repaired and stiffness rises steadily to point E. 

The opposite will occur at the diagram portions marked EF, FG, and GH. After passing the 
minimum angle-of-twist of the unloading portion (EF), progressive resolution of the damaged 
structural components will occur, leading to a slight, unsteady increase in stiffness as negative torque 
increases (FG). After point G, stiffness will increase steadily to point H, as negative torque 
increases. In particular, pinching would correspond to segment CD on the negative torque side, 
where repair of the lesions occurs, and to segment FG on the positive torque side, where resolution 
of the lesions occurs. This explanation, in which lesions form on the negative torque half-cycle and 
reinforcements yield on the positive torque half-cycle, does not take into account buckling. If 
buckling occurs, the situation is reversed: lesions form on the positive torque half-cycle and 
reinforcements yield on the negative torque half-cycle. 

If pinching does not exist, there is no contraction along the cycle. In this case, the slope of 
segment DC on the negative (FG on the positive, respectively) torque side of the model lies in 
between the values of the slopes of ED and CB (EF and GH, respectively), as depicted in Fig. 1 2b. 
The slope of the segments modeling the cycle shows that stiffness decreases along all of them; which 
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might mean that lesions do not repair and do not resolve as they would with pinching. In this 
scenario it is reasonable to assume that torsional loading creates lesions distributed in a more 
disorderly fashion in the osteon than in tensile-compressive loading (aligned with longitudinal 
bundles and where pinching is present). Before load reversal, the hydroxyapatite crystallites might 
have detached and cracked in a way that the original alignment is destroyed. In this case a partial 
realignment does not occur later along the first half-cycle as to bring an increase in stiffness. 
Consequently, there are no lesions to resolve during the successive half-cycle and stiffness would 
keep decreasing. This explanation, in which lesions form on the negative torque half-cycle and 
reinforcements yield on the positive torque half-cycle, does not take into account buckling. If 
buckling occurs, the situation is reversed. Before load reversal, collagen bundles might yield in a 
way that loses track of the original alignment. Similarly this does not allow even a partial 
realignment to take place later along the first half-cycle, after lesions involving the cracking of 
hydroxyapatite have occurred, so as to increase stiffness. Consequently, no partial alignment of 
bundles could be restored before bundles start yielding during the successive half-cycle and stiffness 
would keep decreasing. Such interpretation, where lack of hydroxyapatite crystallite alignment is 
greater under torsional cyclic loading than tension-compression cyclic loading, could be verified 
additionally by means of X-ray diffraction (Ascenzi A. et al, 1 998). 

The results of the measurements and plots above will give the positions and inclinations of the 
segments in Fig. 12. For instance, the t-test on the means of k m along the first half-cycle for 
longitudinal and alternate osteons is expected to indicate that k m is significantly higher for alternate 
osteon. This means that point U in Fig. 12a is higher for alternate osteons, that is, cracks close and 
stiffness switches start increasing for a smaller torque value. This is consistent with the more 
complex structure of alternate osteons. 

D. Fracture Model 

/The purpose of the fracture model is to show that cumulative micro-cracking, de-bonding, 
void growth and fiber breakage associated with repeated loading of osteons causes a progressive loss 
of stiffness and pinching, and increase of energy absorption. The lesions observed under an optical 
microscope in osteon samples subjected to cyclic torsional loading will serve to develop osteon 
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models and to formulate biological hypotheses on propagation of fractures. The fracture model will 
be based on: 

• hypotheses on ultrastructural components' behavior under cyclic torsional loading formulated 
from the experimental hysteretic plots; 

• ultrastructural components' mechanical properties; 

• fractures observed in osteons during monotonic torsional loading; and 

• fractures observed in macroscopic bone specimens. 

This aspect of the model is an adaptation and extension of the approach of Gupta and 
Bergstrom (1998). The fracture propagation model is a micromechanical bone model that allows 
prediction of the progressive growth of faulting zones, by considering the increased stress 
: experienced in the vicinity of an already highly cracked region. The nucleation of initial damage is 
determined by the assessment of the points more susceptible to fracture. The progressive growth of 
the fault nucleus is considered in a statistical manner by the use of stress enhancement factors, which 
address the increased probability of failure in the vicinity of regions that are already cracked. 

The geometric model of each of the longitudinal and alternate osteon samples before 
mechanical testing consists of a hollow cylinder with coaxial lateral surfaces. Its internal diameter, 
external diameter, and height equal 40 urn, 210 urn, and 500 urn, respectively. Each such hollow 
cylinder presents voids, and about 20% of each such hollow cylinder consists of voids (Piekarski, 
1970) which model canaliculae and lacunae. 

The material model of each of the longitudinal and alternate osteon samples before 
mechanical testing consists of a laminate whose length, width, and height correspond to cylindrical 
shell circumference, thickness, and height, respectively (Fig. 9a). The layers are unidirectional fiber- 
reinforced laminae (Fig. 9b) of the same matrix and fibers. The matrix and fibers are each treated as 
homogeneous and isotropic* The fibers are assumed to be circular in cross-section with a diameter of 
800 A, randomly distributed in the transverse plane and perfectly embedded in matrix. The lamina 
with fiber inclination y is named y-lamina. The elastic properties of matrix and fibers model the 
elastic properties of hydroxyapatite (Katz and Ukraincik, 1971) and collagen (Currey, 1969). The 
matrix occupies 40% of the lamina volume without voids (Bonfield and Li, 1 967). The matrix (fiber, 
respectively) volume decreases (increases, respectively) slightly from inner to outer lamina (Amprino 
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and Engstrom , 1 952; Rho et al , ] 999). 

The longitudinal osteon model consists of 9 longitudinal lamellae of the same thickness. 
Longitudinal lamellae are modeled by alternating +82-laminae and -82-Iaminae (Frasca et al, 
1977). The alternate osteon model consists of 7 transverse lamellae with 5 longitudinal lamellae 
5 layered between them (Giraude-Gille, 1988). The transverse lamella is modeled by the laminar 
sequence [-61.5, -41, -20.5, 0, 20.5, 41, 61.5] (Ascenzi M.-G., 1999b). This sequence is 
subjected to prestress as described in Ascenzi M.-G., 1998a and 1999b. A longitudinal lamellar 
model is 9.45 \xxn thick and a transverse lamellar model is 5.40 (im thick (e.g. Gebhardt 1906; 
Ziv et al, 1996). The matrix volume is 1 0% higher in the longitudinal lamellar model than it is 
10 in the the transverse lamellar model (Marotti et al, 1994). 

E. Bone Structure Simulation 

To model fracture propagation in osteons, each of the longitudinal and alternate osteon 
models is divided into a discrete number of elements, e.g. 618,317. The element mesh will be 
15 refined to achieve a convergent solution. A computer simulation program, such as a Monte Carlo 
simulation, will be used to do the following tasks. 

1 . For any given value of torque applied to experimental samples, the distribution of stress in the 
osteon model is computed. Such computation will take into account voids. 

2. Such stress distribution will be added to the distribution of prestress. 

20 3 . The strain associated to the resulting stress will be computed on each phase within each element. 

4. From the strain associated to the resulting stress, the overall deformation of the hollow 
cylindrical shell will be computed. 

5. From the strain in each phase within each element, the phase deformation will be computed. 

6. The strain in each phase within each element is compared to the yield strain. 

25 7. The strain is chosen as the criterion for osteon failure (Piekarski, 1 970). The maximum strain, 
called critical strain, after which fracture occurs within each phase is provided by the literature. 
Perfect bonding at the interface between phases is assumed unless experimentally observed 
cr^si^pcar to initiate at this interface. If that is the case, a failure criterion (e.g. Von Mises) 
will be included. 
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8. The elastic properties of fractured phases will be computed by means of formulas of the type 
ErE c /(l+(Hv c )(k c X)/2) (Gupta and Bergstrom, 1998). 

9. The element is declared broken i F all phases in that element have failed. 

10. The elements are assumed to be aligned in independent rows such that the problem of fracture 
5 propagation becomes one dimensional. 

1 L To model the progressive growth of damage, the torque will be increased incrementally, and 

using the fracture criterion above, the number of failure elements will be established. 
12. The increased probability of fracture in the neighborhood of an already fractured element will be 

considered using the concept of stress enhancement factors. 
10 13. If all elements on one row are broken, the strain level at which all elements on one row are 

broken is taken as the failure strain for that row. The process is repeated for each row in the 

model. Once the maximum torque is reached, the program stops. 
14. At this point, the simulation of fracture propagation before the first hysteretic half-cycle is 

completed. 

15 15. The program incrementally applies a clockwise torque decreasingly to the maximum torque 
applied experimentally to samples, and at each increment it repeats steps 2 and 3 above, so as to 
complete the failure simulation during the first half-cycle of the first cycle. 
1 6. Step 1 5 is repeated for the corresponding counterclockwise torque so as to complete the fracture 
simulation during the second half cycle of the first cycle. 

20 17. The fractures obtained according to the model should be comparable with those observed in 
osteon samples submitted to one cycle only. 

1 8. The fracture simulation sequence is repeated through the last cycle. 

19. Fracture simulation is repeated further, as desired. Fractures obtained in this way, according to 
the model, should be compatible with those observed in cycled osteon samples at the last cycle. 

25 ■ 

F, Results 

Entities computed from the experimental hysteretic diagrams, such as stiffness degradation, 
pinching degradation and increase in energy absorption, up to the second cycle, correlate with the 
fracture propagation of the fracture mode!. 
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The same geometric/material models and computer program will be used to simulate fracture 
propagation under tension, compression, and shear, separately. The resulting fractures will match the 
fractures observed in osteon propagation under tension, compression, and shear respectively 
(Ascenzi A. and Bonucci, 1967 and 1968; Ascenzi A. et aL, 1972). 
5 Predictions and phenomena simulated according to' the model of the invention include that in 

both osteon types, a fracture starts at a weaker point of the bone structure (Carter el ai s 1 981), at the 
weak interfaces between two outer lamellae (e.g. Piekarski, 1970), presumably because of the 
hydroxyapatite decrease in osteons from vascular canal to outer wall (Rho et aL, 1999). 

In longitudinal osteons, the fracture starts somewhat longitudinally, between collagen 
10 bundles. It then deviates once or twice at the beginning of the fracture and is soon followed by a 
smooth crack advancing rapidly across the osteon to possibly end in the vascular canal. As cycling 
continues, collagen bundles between cracks break, and cracks join to create one or more long almost 
vertical cracks. 

In alternate osteons, cracks are expected to spread obliquely by following the weak interfaces 
15 of lamellae. The transverse and oblique collagen bundles may break before the longitudinal ones as 
the osteon section enlarges. Cracks spread through lamellae less rapidly than in longitudinal osteons 
as explained by the crack propagation control, characteristic of composite materials (Cook and 
Gordon, 1964). Once the crack breaks through transverse and oblique bundles, it will propagate 
faster straight through the vascular canal. A long crack should show an oblique orientation between 
20 upper and lower extremities. 

It follows, unexpectedly, that the longitudinal osteon is weaker in longitudinal than in 
transverse shearing, while expectedly the alternate osteon is weaker in tension than in shear (Ascenzi 
et aL, 1 967 and 1 972). This is because when a torque is applied to a body, tensile and compressive 
stresses are produced on the lateral surface and torsional shearing stresses are produced on the cross- 
25 section of the body. The tensile and compressive stresses act approximately at a 45° angle to the 
long axis of the body. The direction of the shearing stress on the cross-section of the body is the 
same as that of the force producing torsion. If a material is weaker in longitudinal than in transverse 
shearing, the first cracks arise from axial shearing stresses and appear in a longitudinal direction. 
However, if the material is weaker in tension than in shear, it usually cracks along a spiral course 
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inclined at a 45° angle to the long axis of the body. The reason for this is that a state of pure shear is 
equivalent to a state of tension in one direction and of compression in the opposite direction 
(Timoshenko and Young, 1962). The tension stress produces a spiral crack in the body. 

For both osteon types 3 to 4 small cracks form in the hydroxyapatite and collagen, which 
5 yields and pulls and/or buckles and makes the cracks spread within lamellae. Microcracks form 
ahead of the advancing fracture line. Afterwards, during torque reversal, width of cracks and 
deformation decrease. Collagen may buckle and some resisting strength may appear at zero torque. 
As cycling continues, cracks extend through the lamellae and join. 

The slow propagation of cracks in the areas containing transverse and oblique collagen 
10 bundles allows for the area to absorb a large amount of energy. Slow propagation is essentially a 
pull-out type mechanism, that is, hydroxyapatite crystallites are pulled out of the collagen by shear 
failure at the fiber-matrix interface. The rapid propagation of cracks in areas containing 
approximately vertical collagen bundles allows very low energy absorption. This should be 
compatible with larger areas enclosed by cycles of experimental plots of alternate osteons (see last t- 
15 test in Sec, A above). 

Hydroxyapatite crystallites are pulled out from collagen around canalicular 

At low strain rates in compression distortion of the lamellar structures occurs (McElhaney 
and Byars, 1965). 

The propagating crack generally has the tendency to avoid discontinuities (Piekarsly, 1 970) 5 
20 hence increasing its length. Discontinuities act as crack arresters by blunting the tip of the crack 
which enters them. 

;v The fracture model is expected to agree with fractures observed in osteons cycled only for 

first and second complete cycle. The dimensions of the hollow cylindrical model after one cycle, 
two cycles, and the last cycle of torsional loading should match the means of the osteon samples 1 
25 dimensions measured experimentally. Furthermore, the role of the models' fibers is expected to 
check with the cyclic behavior of decalcified osteons. 

The sudden shift of the osteon shape (Fig. 4) from a circular to a square cross-section 
suggests a stress concentration at the lugs. Therefore, fractures may begin at the end of some 
samples earlier during loading than would otherwise be expected. 
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Lamellar thickness and width were measured on 20 bright and 20 extinct peripheral lamellar 
samples in quintuplicate in dry osteon samples by Delta Sistemi IAS 2000 image analysis system, 
and again after wetting with a micro-pipette. This table shows means and standard deviations. 
Thinner extinct lamellae were used for comparison with bright lamellae. It is known that extinct 
5 lamellae are thicker than bright ones, whether dry or wet. The Student t-test is run on the data to 
determine statistical differences between dry and wet lamellar dimensions. 



Sample 


Thickness 


Thickness 


Width 


Width 




Dry 


Wet 


dry 


wet 


Bright 


3.30+0.88 


3.56+0.93 


70.30±9.28 


72.45+9.58 


Extinct 


4.13±1.23 


4.10±1.10 


70.30±9.28 


72.45±9.58 



Whether dry or wet, bright lamellae are significantly thinner than extinct lamellae when 
enclosed in alternate osteons (this agrees with previous results, e.g. , Rho et al. s 1 999). Additionally, 
wet and dry conditions affect bright and extinct lamella thickness differently. Bright lamellae are 
1 0 significantly less thick when diy than wet. In contrast, extinct lamellae thickness remains constant 
whether wet or dry. 

The bright lamella thickness increases from dry to wet which may be due to the higher quota 
of mucopolysaccarides, which expand with water, and to the transverse collagen bundles in the 
bright lamella tightly encircling extinct lamella, thereby impeding expansion. Height of both bright 
1 5 and extinct lamellae is significantly smaller when dry. In addition, the thickness along lower and 
upper borders shows variations up to 50-60%. This will be included in the model. Width variation is 
very low. 

The model provides an advantageously simplified simulation or representation of osteon 
structure. For example, partially calcified collagen bundles are excluded from the model. The model 
20 provides a useful and improved description of bone structure and mechanics, even though the shape 
and dimensions of hydroxyapatite crystallites and the relationship of these parameters to the organic 
components of the matrix are only partially known. Not all collagen bundles are completely 
calcified. Those, which are note calcified take up crystallites only on 400 A bands (Ascenzi, A. et 
al, 1 965). Such bundles may be comprised of relative more stiff 400 A bands separated by relatively 
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more flexible decalcified collagen segments. In a preferred embodiment of the invention, partially 
calcified collagen bundles are not modeled, in favor of modeling fibers in uncalcified collagen 
bundles. The matrix, which models the hydroxyapatite crystals, lies outside the fibers. Pinching is 
incorporated into the model is related to the yielding and bucking of fibers, and provides an 
5 approximation of the yielding and buckling of partially calcified collagen bundles. In preferred 
embodiments, fracture propagation is modeled and cracks will tend to propagate before buckling is 
likely to occur, because the model in most cases assumes that individual fibers are perfectly bonded 
to and are uniformly supported by the matrix. The model also excludes complex consideration of 
pore fluids in preferred embodiments which balance relative simplicity with achieving a reliable and 
10 accurate bone model. 



The patents, applications, test methods, and publications mentioned herein are hereby 
1 5 incorporated by reference in their entirety. 

Many variations of the present invention will suggest themselves to those skilled in the 
art in light of the above detailed description. All such obvious variations are within the full 
intended scope of the appended claims. 
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Wc Claim : 

1 . A model of macrostruciural properties of a bone, wherein said model comprises 
hierarchical structural and hierarchical mechanical properties of microstructure of said bone and 
interactions of said bone with external force. 

2. A model as defined in claim 1 , wherein said bone is compact bone or cancellous 

bone. 

3 . A model as defined in claim 1 , wherein said mechanical properties are selected 
from the group consisting of tension, compression, shear, bending, torsion, prestress, pinching, 
and cement line slippage. 

4. A method of predicting deformation and fractures of bone using the model as 
defined in claim 1 . 

5 . A method of identifying the requirements of bone reconstruction and prosthesis 
using the model as defined in claim 1 . 
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